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Using Galaxy Cluster Peculiar Velocities to Constrain Cosmological 
Parameters 

A. Peel a * and L. Knox at 

a UC Davis Cosmology Group, 

1 Shields Ave, Davis, CA 95616, USA 

Galaxy cluster peculiar velocities can be inferred from high-sensitivity, high-resolution multiple-frequency 
observations in the 30 to 400 GHz range. While galaxy cluster counts and power spectra are sensitive to the growth 
factor, peculiar velocities are sensitive to the time-derivative of the growth factor and are hence complementary. 
Using linear perturbation theory, we forecast constraints on fl m , Ho and the dark-energy equation of state 
parameter, w, given 820 densely sampled cluster locations (at z ~ 1) from a ACDM N-body simulation and 820 
sparsely sampled cluster locations in a broader redshift range. 
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1. Introduction 

After the cosmic microwave background, 
galaxy clusters may be our best hope for a 
precision cosmological probe. Although clus- 
ters are highly non-linear objects, the domi- 
nant importance of gravitational effects over non- 
gravitational ones makes them attractive can- 
didates as systems whose properties should be 
calculable from first principles. Cluster number 
density evolution has been proposed as a way 
to determine cosmological parameters including 
the dark-energy equation of state parameter [Q . 
Here we investigate constraints on cosmology pos- 
sible from measurements of cluster peculiar veloc- 
ities. 

Peculiar velocities of galaxies measured via, 
e.g., the Tully-Fisher and fundamental plane re- 
lations have already been used to constrain cos- 
mological parameters Q . These efforts have been 
hindered by large statistical and systematic errors 
in the velocity measurements. In contrast, clus- 
ter peculiar velocities determined from the kinetic 
and thermal Sunyaev-ZePdovich (SZ) effects are 
potentially more reliable. In addition, SZ is the 
only way to determine peculiar velocities over a 
wide range of redshifts allowing one to observe 
the evolution of statistical properties. 
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2. Sunyaev-Zel'dovich Effects 

About 1% of the cosmic microwave background 
(CMB) photons that pass through the gravita- 
tionally confined hot gas in the deep potential 
wells of galaxy clusters Compton scatter off a free 
electron. On average these photons gain energy 
causing a spectral distortion of the blackbody na- 
ture of the CMB characterized by a deficit of low 
frequency photons and an excess at high frequen- 
cies, with a null near v = 217 GHz. Specifically 
the distortion has the frequency dependence || : 
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This spectral distortion has been observed in the 
direction of numerous known clusters, whose lo- 
cations had been determined by other means (op- 
tical and/or X-ray surveys) Q. 

The bulk motion of the cluster with respect 
to the CMB causes the kinetic SZ effect which, 
to first order, leads to a hotter (colder) black- 
body spectrum for a cluster that is moving toward 
(away) from the observer ||: 
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where v r is the radial component of the velocity. 
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Measurements at multiple frequency bands in 
the 30 to 400 GHz range can be used to simulta- 
neously determine y ~ r J^^-i and Tv r /c. These 
two determinations together with a temperature 
determination can be used to solve for v r /c j|: 
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For hotter clusters there are temperature- 
dependent corrections to Eq. [l] that allow one 
to solve for the temperature as well. For cooler 
clusters, X-ray determinations of the tempera- 
ture may be necessary. 

The best determinations of peculiar velocities 
from SZ measurements have errors of a v ~ 1000 
km/s H but can be greatly reduced with more 
sensitive measurements. Published predictions 
of achievable a v [|J assume that the motions of 
the electrons in the cluster are simply thermal 
motions plus one coherent bulk motion. In this 
case there is indeed a well-defined "velocity of the 
cluster", and the error in its measurement is: 
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In the no-noise (AT noisc = 0) limit and with a 
confusion noise from the Ostriker-Vishniac effect 
of ATcmb = 2 /iK (achievable for clusters of 
small angular extent which includes any cluster at 
z > 0.2, and with high (~ 1') angular resolution) 
then a v ~ 25 km/sec is possible. 

However, clusters have many bulk flows in 
them, not just one. Integrating over all these in- 
ternal bulk flows, with the mass-weightings natu- 
rally provided by the SZ effects, should still lead 
to what might be called "the velocity of the clus- 
ter" . But a sizeable fraction of the mass is along 
lines of sight with quite low optical depth and 
hence high velocity errors. The net result is that 
a v ~ 100 km/sec is a more likely error M. 

3. Survey Strategy Issues 

We imagine peculiar velocities determined as 
part of a high-resolution, multi-frequency (30 to 
400 GHz) follow-up campaign on a cluster survey. 
These clusters may have been found from X-ray, 
SZ or even galaxy redshift surveys (e.g., SDSS 



and DEEPII). Important issues are how and how 
accurately to determine redshifts, how many ve- 
locity measurements to make at each redshift and 
whether to sample densely to study correlations 
or sparsely to cover more volume. 

Cluster redshifts may be obtained via (i) op- 
tical photometric and/or spectroscopic redshifts 
of galaxy cluster members, (ii) measurement of 
CO lines in galaxy cluster members, or (iii) X- 
ray measurements with sufficient spectral resolu- 
tion. X-ray and optical data also provide valuable 
information on the dynamic state of the cluster; 
e.g, has the cluster fully virialized, or is it still 
merging? Cluster redshifts are already vitally im- 
portant for determining dN/dz from thermal SZ 
surveys. Note that interpretation of peculiar ve- 
locity correlations is much less sensitive to selec- 
tion effects than the dN/dz /dQ, statistic. 

Dense sampling puts much higher demand 
on the accuracy of the redshift determinations. 
For a specified error in comoving separation Ar 
we need a redshift error smaller than Az = 
0.003/i(#(z)/i7 )(Ar/(10 Mpc)); in order to cal- 
culate the expected correlations accurately we 
need Ar < 10 Mpc. With sparse sampling, all 
we need is a coarse redshift binning and assur- 
ance that the clusters are separated by > 200 Mpc 
(comoving); photometric redshifts are almost cer- 
tainly sufficient for that case. 

There are good arguments for targeting the 
z ~ 0.2 to 0.5 range rather than higher red- 
shifts. The tolerance on redshift error decreases 
only slowly with increasing redshift, whereas the 
amount of telescope time required typically in- 
creases dramatically. In addition, optical and X- 
ray data on the clusters may be valuable for in- 
vestigating the dynamical state of the cluster and 
third X-ray temperatures are necessary for the 
colder clusters. 



4. Linear Theory 

The evolution of the Fourier-transformed den- 
sity contrast 5k in the linear regime is separable, 
allowing us to write: 



5k(ri) = D(r))5 k (r)a) 



(5) 
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where D(rj) is the growth factor and ijo is the con- 
formal time today. From the continuity equation 
for matter, ikvk = Sk, we see that velocities are 
probes of the time derivative of D. 

The correlation between radial velocity com- 
ponents of two clusters at locations (rj,7i) and 
(rj , 7j ) relative to the observer is given by || : 
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where cos# = % ■ jj, K±(kr) = ji(kr)/(kr), 
K\\(kr) = ja{kr) — 2ji(kr)/(kr), r is the comov- 
ing distance between the clusters and an overdot 
symbolizes d/drj. Figure Q shows ^f±(6) for four 
different flat cosmologies at redshift z = 1. 

Cosmology dependence in Eq. |^ arises from: 
(i) the redshift-distance relation r{z) and (ii) the 
time-derivative of the growth factor, D. The lat- 
ter has weak dependence on w as we can see from 
the highly accurate analytic approximation JfTJ]: 



D{z) = D(z)a(z)H(z)[n m (z)] a ^ 



(7) 



because a only ranges from about 6/11 to 3/5 for 
w = — 1 to w = 0, and because a(z)D(z)H{z) is 
also fairly insensitive to w for z < 0.5. In con- 
trast, D is highly sensitive to f2 m . 

5. Forecasted Parameter Errors 

For nearly uncorrelated cluster velocities, from, 
e.g., a very sparse survey, we can easily estimate 
the expected error variance on f2 m from measure- 
ment of N clusters with expected peculiar velocity 
variance ^o{zi): 
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where the last equality assumes N clusters with 
cr 2 < all at z = 1, and d\n^ /dfl m ~ 5. 
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Figure 1. Linear peculiar velocity correlations for 
pairs of clusters at redshift z = 1 MS M. function of 
separation angle. 



We now turn to a Fisher matrix analysis of 
two survey types. For a dense survey, we sample 
820 neighboring clusters from the Hubble Volume 
Lightcone Simulation Cluster Catalog iJTTJ] over an 
area of 80 sq. deg. at redshift z ~ 1. For a sparse 
survey, we construct a grid of 820 clusters from 
z — 0.1 to 1.0 over an area of 2600 sq. deg such 
that all comoving cluster separations are greater 
than 200 Mpc. 

We calculate the expected covariances in linear 
theory for a density field smoothed with a 5 Mpc 
radius tophat filter. We calculate their partial 
derivatives by finite difference. The models we 
difference have fixed COBE normalization p2[ . 
All our results in Figure ^ assume a noise level of 
100 km/s. 

As expected, fl m is measured very well. Nei- 
ther h nor w are well constrained; nor do they 
lead to significant confusion in the f2 m measure- 
ment. In addition the results are very similar for 
our two different survey types, with the sparse 
survey (dashed lines) marginally more constrain- 
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Figure 2. Forecasted 1 and 2-cr contours for the 
dense survey (solid) and the sparse survey with 
(dN/dz)Az as shown (dashed). For each two- 
parameter contour, the third parameter is held 
fixed. [The drop in dN/dz in the final bin is due 
to a quirk in our sampling algorithm.] 



ing than the dense survey. Both increasing the 
velocity errors to greater than the assumed 100 
km/sec and including the effects of redshift er- 
rors will improve the sparse sampling result in 
comparison to the dense sampling result. 

6. Non linear and biasing effects 

We expect that linear theory provides only 
a rough guide to the constraints possible from 
peculiar velocity measurements. Numerically- 
determined variances agree with linear theory 
predictions to ~ 40% in the present epoch |X3| . 
We plan to include non-linear and biasing effects 
by use of the P1NOCCHIO gravitational instabil- 
ity code Q- With P1NOCCHIO we will obtain 
highly accurate calculations of the velocity corre- 
lation functions in several models and once again 
use finite difference to calculate Fisher matrices. 



7. Conclusions 

Cluster peculiar velocity measurements are 
highly sensitive to fl m and largely insensitive to h 
and w. Sparse sampling is probably preferable to 
dense sampling. Velocity determinations of ~ 800 
clusters can potentially be used to determine fl m 
to ~ 0.03%. This w-independent measurement 
of Q m may be useful for improving w constraints 
from observations of type la supernovae |l5[ ] . 

We would like to thank B. Holdcn, G. Holder, 
M. Kaplinghat, C.-P. Ma, J. Mohr, S. Meyer, M. 
White, G. Wilson and I. Zchavi for very useful 
conversations. 
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